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PRODUCTS OF MATRICES 


1 1 
0 1 


AND 


1 0 
1 1 


AND THE 


DISTRIBUTION OF REDUCED QUADRATIC IRRATIONALS 


FLORIN P. BOCA 


Abstract. Let $(IV) denote the number of products of matrices [J (] and [( °] of trace 
equal to N, and 'F(iV) = Yl n =3 ^( n ) be the number of such products of trace at most N. 

We prove an asymptotic formula of type 'fr(N) = ciN 2 log N + C 2 N 2 + O e (N‘ / 4+e ) as 
N —* oo. As a result, the Dirichlet series ^(n)n~ s has a meromorphic extension 

in the half-plane 3t(s) > | with a single, order two pole at s = 2. Our estimate also 
improves on an asymptotic result of Faivre concerning the distribution of reduced quadratic 
irrationals, providing an explicit upper bound for the error term. 


1. Introduction 


Consider the matrices 


A = 


1 

1 


0 

1 


and 



1 

1 


and the (free) multiplicative monoid Ai they generate. The problem of estimating 


$(N) = #{C e M : Tr(C) = N} 


q>(N) = $(n) = #{C G M : Tr(C) < N} 

n =3 

for large N came across in the study of a number-theoretic spin chain model in statistical 
mechanics introduced in ms, and further investigated in [5] and [2j. In [2j the estimate 

N 2 log N o 

(1.1) *(N)= c(2 ° + 0(N 2 log log N) (IV—► oo) 

was proved, using L -functions and a result from ,1J concerning the distribution of reduced 
quadratic irrationals. 

In this paper we improve the estimate a. Our approach relies on a result concern¬ 
ing the distribution of multiplicative inverses, which is a consequence of Weil’s bound on 
Kloosterman sums. 
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Theorem 1.1. We have 

(L2) — 
where 

( 1 . 3 ) MN) < £ Nl +£ 

and 


V(N) = ci-ZV 2 log N + c 2 N 2 + $ 0 (JV), 


(iV -> oo), 


1 


1 


Cl = 


C2 = 


^ib l0E2 -ff} 


C(2) ’ ^ C(2) 

Using the Mellin transform representation of Dirichlet series and changing N to N — 2 
in the right-hand side of m we obtain for 3i(s) > 2 




\ — s —1 


= S 


POO 

/ \k(x)(x — 2) 

./3 

/ OO POO 

{c\X 2 logX + C 2 X 2 )x ~ S ~ 1 dx + S / 

OO 


dx 


2ci 

(s-2) 


, C! + 2 C 2 
H---1- C 2 + s 


s -2 


Since the function \ko(aO satisfies the growth condition D the integral above converges 
for 3i(s) > |, and thus the right-hand side defines an analytic continuation of Z{s) to the 
half-plane 3?(s) > ^ with the point s = 2 removed. 

The contribution to the main term in the asymptotic formula above only comes from 
words of odd length in Xi. We estimate the contribution ^f ev (N) to 4 f(N) of words of even 
length which begin in A and end in B , proving 


Proposition 1.2. 


^ ev(iV) = ^(^ + 0£(iVl+£) 


(N 


oo 


This extends Proposition 4.5 in P| which states that 

T /A 7 \ IV 2 log 2 
ev( - ^ ~ 2C(2) ’ 

without any estimate on the error term. 

Using a transfer operator associated with the Gauss map, Fredholm theory, and Ikehara’s 
tauberian theorem, Faivre |lj proved the asymptotic formula 

e x log 2 


E i 

p(u>)<X 


2C(2) 


(X 


oo 


for the number of reduced quadratic irrationals lo of length p(ui) at most X. Since the final 
argument relies on a tauberian theorem, no explicit bound was found for the error term. 
In the last section we use Proposition 11.21 and an explicit identification from between 
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products of matrices A and B starting with B and ending with A, and reduced quadratic 
irrationals, to prove 

e A log 2 


Proposition 1.3. 1 

p(u>)<X 


2C(2) 


+ O e (e ( 8+ £ ) x ) [X —► oo). 


If 


2. Products of ^4’s and B' s and continued fractions 


a — [ai, a 2 ,... ] — 


CL 1 + 


a 2H--T- 

“3+t!t 


( 2 . 1 ) 


is a reduced continued fraction with positive integers cu, the k m convergent 

— = [ai,... ,a k \ 

Qk 

is given by pairs (PmQn) of relatively prime integers defined recursively as 

PO = 0 , Pi = 1 , Pn = a n p n - 1 +Pn- 2 , 

qo = lj 9 l = ^ 1 ? Qn — ^nQn —1 H” QVi— 2 ? 

and satisfying 0 < p n < q n and the equality 

Pn—lQu PnQn—l = ( 1 ) • 

If (p n ) and (g n ) satisfy (|2.1() for every 0 < n < k, then 

qn— 1 qn— 1 1 


QVi ttnqn—1 H” qn—2 &n H - 


Qn — 2 ’ 

qn-l 


showing that 

( 2 . 2 ) 

Consider the matrices 


gfc-i 

<7fc 


, • • • , Ul] • 


J = 


0 1 
1 0 


and M(a) = 


a 1 
1 0 


qn qn —1 
Pn Pn —1 


If ( p n , gvi) is as in (12.11) then, as noticed in JJJ, we have 

(2.3) M{a 1 )... M{a n ) = 

When combined with 

B k A e = k,i€ Z, 

equality yields 

(2.4) B ai A a2 • • • B a2m - x A a2m = Af(ai) • • • M(a 2m ) = 


(Z27?i Q2m—1 

P2m P2m—1 
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From EH, EH), and 


we also infer that 


B = A T = JAJ , A = JBJ, 


A ai B a2 • 

. J^2m-1 g«2m _ J 

Q2m 

<72771-1 

J = 

P2m — 1 

P2m 



P2m 

P2m — 1 


92m-1 

92771 


(2.5) 


B ai A a2 ■ ■ ■ A a2rn B a2m+1 = 


92m 92m-1 

P2m P2m—1 


1 «2m+l 

0 1 



0.2m 

92m+l 


P2m 

7*2777+1 


jyll 1^0,2 _ _ _ J^a2m J^a,2m+l _ 


7*2/77+1 7*2771 

92m+l 92m 

All matrices in the products from EH and (12.51) have determinant 1. We denote 
W ev (N) = {(ai,...,a 2m ) G N 2m : m > 1, Tr(£> ai A a2 • • • B a2rn ~ 1 A a2rn ) < N}, 
W od d(A0 = {(ai,... , 02 ^+ 1 ) G N 2m+1 : m > 1, Tr(B ai A a2 •. < AT}. 

We consider the sets 


S ev (N) = 


q q 

p' p 


0 < p < q, 0 < p' < q', q' > q, 
p + q' < N, pq' — p'q = 1 


and 


Sodd(iV) = 


q q 


0 < p < q, 0 < pf < q', q' > q, 1 
p p'j p' + q < N, p'q — pq' = 1 J ’ 

of cardinality fy ev (N) and respectively ^ 0 dd(-^), and the maps defined as 

/3ev(®l, • • • , ®2 m) — M(ai) • ■ ■ Mi y CL2m)i /3odd(®l, • • • , ®2m+l) — A/(oq) • • • Af((l2m+l): 

from Uf 3 W ev (A') to U ^°S ev (N), and respectively from U^Woddl-^V) to Uf 3 5 0 dd(A r ). As 
a consequence of EH, /3 ev and /3 odd are injective. It follows from EH and (12.51) that 
/3ev(W ev (A0) c S ev (N) and P odd (W odd (N)) C S odd (^)- To check S ev (N ) C /3 ev ( W ev ( N) ), 


let 


<? Q 

v' v 


G S ev (N). With K = [32] we have 0 < q' — Kq < q, 0 < p’ — Kp, and 

(j)' — Kp)q— (q' — Kq)p = 1. Since q' > q are relatively prime, we also have K < , 

and thus p' — Kp < q' — Kq. Since 


q q 

p' p 


0 1 
1 —K 


q q 

p' p 


M(K)- 1 = 


q q' — Kq 
p p' — Kp 


it follows, by replacing (q. q’) by {q' — Kq , q) and (p,p') by (p 1 — Kp,p) and performing this 
process until q' — Kq becomes equal to 1, that the matrix q is written as a product 

of k matrices of form M(K). Since q' > q, k ought to be even and therefore p , p G 
/3 ev (W ev (N)). One shows in a similar way that S odd ( N) C /3 o dd0 / V od d(AQ). 
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This proves that the elements of M. are uniquely represented as products of A’s and B' s. 
It also implies that 

(2.6) T(1V) = 2\H ev (lV) + 2* odd (N). 


3. Estimating T ev (A0 


To estimate T ev (A^), we first keep q' and q constant. From pq' — p'q = 1 and q < q', 
it follows that q and q' are relatively prime, and that p is uniquely determined as p = q', 
where q' is the unique integer in {1,..., q} for which q'q’ = 1 (mod q). It is obvious that 
p’ := pq ~ 1 < q' and the map 


(q,q') ■ 


q < q’ < IV, (q, q’) = 1 
q' + q' < N 


s (■ q,q') 


q 

j _ „ _ 


p = 


p = q' 


G S ev (N) 


is a one-to-one correspondence. Replacing q' by y and q' by x we can write 


(3.1) *ev(A0= E 1 = E E L 

q<q'<N q<N q<y<N 

qf+q’^N 0<x<min{q,N-y} 

xy= 1 (mod q) 

For each y E (0, N], there is at most one x E (0, q) such that xy = 1 (mod q); whence 
the trivial estimate 

T ev (IV) < N 2 . 

To give a more precise estimate for \E , ev (IV), we shall define for q > 1 integer and subset 
in M 2 the number of (relatively prime) integers (x,y) E II such that xy = 1 (mod q). 

It is known (see for instance |T]) that Weil’s bound on Kloosterman sums yields for any 
intervals I and J of length at most q the estimate 

AT q (IX j) = ^f \I\\J\+O e ( q ^). 

qz 

This immediately extends to intervals of arbitrary size as 

(3.2) A f q (I xJ) = ^ |/| \J\ + 0 £ (^ +e (l + 7 ) ( X + y)) ' 

Another easy but useful consequence of (Em is given next. 


Lemma 3.1. Suppose that 

M = {( x,y ) : a<x<P, fi{x) <y< ^(x)}, 
with C 1 -functions fi < fi on [a,/3]. For every positive integer T we have 

J\fq{Q) = Area(fl) + £ q , 
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with 


£q 


(3 — a 

~w 


{V£{h) + V£{h))+Tq^(l + 


(3 — a 

~w 


1 + 


WflWoo + II /s lie 


where Va{fi) denotes the total variation of fi on [a,(3\. 


Proof. One can take without loss of generality /i = 0 and /2 > 0. Partitioning [a, (3\ 
into T intervals Ii of equal size and denoting by Mi and rrii the maximum, respectively the 
minimum, of fi on /y, we clearly have 

T—l T—l 

X [°> m *]) <-^(0) < X [°> M i ])• 

z=0 z=0 

By (Id.21) we infer that 


Nq{k x [0,mj]) = 


¥>(?) P-a 


rm + O e [ q 2 +e ( 1 + 


/3 — a 

w 


1 + 


11/2 ||c 


The statement follows from this, the similar estimate for M q (li X [0, M*]), and from 


rP T-l 

/ f2(x)dx = Y^ 

Ja 1=0 


(3 — a 


mt + O 


T-l 


^ V’w) = E ^ + o 


i=0 


/3 — a 


Vi( A) . 


□ 


Equality EU can also be written as 
(3.3) 


*ev(A0 = XI -WiV, 

q<N 


qji 


where f^N, q = {(a?, 2/) : q < y < N, 0 < x < mm{q,N — y}} coincides with the trapezoid 
^Tvly = {{ x iV) '■ 0<x<q<y<N — x} if q < and with the triangle = {(x, y) : 
0 < x < N — q, q<y<N — x} if q>Y- Employing (13.31) and applying Lemma ITHI to 
and we infer that 


T ev (iV) = XI Af q (Cl§ q ) + XI ^n] 


)( 2 ) 




-<q<N 


E ^|).(L^)f + £ l(W)+ £ 2(W) , 




%<q<N 



with 
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Taking T = [N 4 ] and using standard summation results (see formulas (14.61) . (14.71) . (14.81) 
below) this gives 

*eAN) = ^^ + o £ ( m + t ), 

which ends the proof of Proposition 11.21 


4. Estimating 'S'odd(N) 


To estimate ^oddi^), we first keep q and p fixed. The general solution of p'q — pq' = 1 
is given by 

/ - , . / P'Q ~ 1 QQ~ 1 . .-r. 

p = q + pt, q = - =-h qt, t E £, 


p 


P 


where q is the unique integer in {1,... ,p} for which qq = 1 (mod p). Since p' > p, one has 
t > 1. The map 


(Q,P,t) 


0 <P<q, ( P,q) = 1] 


3 {q,p,t) 


Q 

p 


q i = 2tzi + qt 
p' = q+pt 


G ‘Sodd (^V) 


is a bijection. Replacing p by a, q by y, and q by x, we can write 


^odd(iv) = £ £ 

q<N p<q 
q+q<N 


N — q — q 
P 


£ 


E 


a<N a<y<N 

0<x<min{a,N—y} 
xy= 1 (mod a) 


N — y — x 
a 


When N — y < a we have N — y — x < a, and the contribution of such terms to 4 , 0 dd(-^' r ) is 
null. So we only consider N — y > a. This gives a < ^ and thus 


*odd(A0 = £ £ 


From 


it follows that 


N — y — x 


a<E. a<y<N—a 
2 0<iK<a 
xy =1 (mod a) 


< 


N — y — x 


N — y — x 


< 


N-a 

N 

< - 




a 

a 

N' 

1 


a 




The set of points (x, y) E (0, a] x (a,N — a) for which [ If y ] = i coincides with 


&N,a,i '■= {(x, y) E (0, a] x (a, N — a) : N — [i + l)a < x + y < N — ia}. 





























We can thus write 


[2]-i 

^odd(iV) = E E iAf a {^N,a,i) 


■n i =1 


(4.1) 


where the sets 


[f]-i 

(j'Mi (^AT,aj) + (j ~ l)-^a(^A/>,:/))’ 


a<f 3=1 


Wv,a,j = {(^ 5 2/) e (0, a] x (N - (j + l)a,N - ja], x + y < N - ja}, 
n( N,a,j = {( x >v) e (0,a] x (IV - (j + l)a, JV - ja], x + y > IV - ja}, 1 < j < 

(i) 


JV' 

a 


- 2 , 




17 


JV,a,[JV/a]-l 
( 2 ) 

N,a,[N/< 


. o <x<N- [^]a, a < y < N - ([^] - l)a, 
x + y<N- ([£] -l)o 


al—1 = i ( X ^) : 


0 < x < a, a < y < N — ([^] — l)a, 


x + y > N - ([f] - 1 )a j ’ 

give a partition of the trapezoid {(x,y) : 0<x<a<y<N — a — x}. Since A/" a (17) does 
not change when 17 is translated by integer multiples of (a, 0), the right-hand side in (14.11) 
can also be expressed as 

E +V a (5 N ,„ ,[JV/a]-l), 

3=1 

where for 1 < j < [^] — 2 we set 

<SN,a,j = {(x,y) G (0, ja] x (iV - (j + l)a,JV - ja], x + y < N - a}, 
and for j = [^-] — 1 we set 


— 1 )a, x + y < N — a > . 


^N,a,[N/a]—i = j0+ y) ■ 0<x<iV-2a, a < y < AT - ^ ^ 

The sets 5jv,aj, 1 < j < [^] — 1, are mutually disjoint and their union is the triangle 
7jv,a = {(+ y) : 0 < x < iV — 2a, a < y < N — a — x}, 

and thus we get 

^odd(JV) = M a {T N , a ). 

a<f 

Since J\f a (17) is invariant under translations by (0, a), this further gives 
(4.2) * odd (A0 = E A f a (X N , a ), 
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where 

TJv i0 = {(x,y) : 0 < x < N — 2a, 0 < y < N — 2a — x} 
is the translated triangle T/v, a — (0, a). 

Lemma 4.1. For every 0 < c < 1 we have 

E = E ^ +o(iv I+c ). 

a<N c a<N c 

Proof. Set K = [A] — 2 > 1. We partition the triangle 7/v,a as T>N,a U ~R-N,ai with 

K 

F>N,a = IJ (o, ia] x ((A' - i)a, (K - i + l)a] 
i =1 

and Ajv, a = T/v,a \ Pn.o,- As T>N,a is the union of 1 + 2 + • • • + K = — disjoint squares 
of size a, we have 

(4.3) A r a (V N ,a) = A(i l +1) <p(a) = ^ AreaVjy a . 

2 a- 

On the other hand it is clear that 

N 

Na{fcN,a) < (K + 1 )ip(a) < — • a = N, 

a 

which gives in turn 


(4.4) 


Y Ma(n N}a ) < N 


1+c 


a<N c 


2 ' a 

a<N c 


We also have 

(4.5) Y AreaA-AT a < Y ~ AreaAW,a < Y --(AT+1)-< Y —-a<N 1+c . 

a<N c a<N c a<N c 

Employing H31), (O and (HU), we gather 

Y M a (T N , a ) = Y M a (V N , a ) + 0(N 1+C ) = Y 


fifl) 


Area V N>a + 0{N 1+C ) 


a<N c 


a<N c 


a<N c 


Y (Area T N>a - Areata) + 0(N 1+C ) 


a<N c 


E 

a<N c 


<p(a) 


• AreaT/v,a + 0(A 1+C ) 


2 ^■ (JV - 2fl)2 + Q(]yi+ c ). 


a<N c 


□ 
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Lemma 4.2. For every 0 < c < 1 and every integer T > 1 we have 

' A 2 log N 




JV c <a<f 


N c <a <f 


^IV 2 log IV + riv f+ £ + Ar 2 +( £ -|)c^ _ 


Proof. Applying TiRnima l.ll I we get 

A fa(f N ,a) = ■ Area (f N , a ) + O e 


N 

Ta 


A + Ta^ +£ ( 1 + 


N\ N 
Ta) a 


<p(a) (A" — 2a) 2 , „ f N 2 


+ O e ( — + TNa~^ +e + A 2 a _ i +e 
Ta 


which is summed in the range a 6 (1V C , y) to get 


V- ¥>(«) (A-2a) 2 , ^ ^ A 2 log A 

2s „2 9 + 


N c <a<% 


N 


— o+ e 


a 2 


E 


¥>(<*) (N-2af . 

o’ ^ i 


N c <a<- 


T +™E« _J+ ' + iv2 E 

+ rjv | +e + N 2 + (.-i)c\ 


□ 


Taking T = [IV 4 ] and (any) c G [ 2 , 3 ], the previous two lemmas, together with (14.21) . 


yield 

Corollary 4.3. ^ odd (A) = ^ 


< d a MJV-2a) 2 +0g ( Jvjn 


“<f 


Lemma 4.4. Sat := ^ 9 -— = Cjy + 0(A), where 


2 a 2 


a< 7 


Otv = 


IV 2 


2C(2) 

Proof. Employing the Dirichlet series 

<p(a) _ C(s - 1) 

CM ’ 


(logiV + 7 -log2-|-£|i) 


E 

a =1 


and the Perron integral formula (< 7 o > 0) 


^ /-o-Q+ioo 


kl 


y 


CTQ —IOO 


s(s + l)(s + 2 ) 


ds = 


s = a + it, 0 > 2 , 


0 if 0 < y < 1 , 

Mi-J ) 2 ify> 1. 


with y = h we infer that 


^ /-cro+ioo 

Sn = / s(«) ds, 


27ri 


' (TQ—IOO 
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with 


9(s) = 


jys+2 


C(s + 1) 


2 s (s + l)(s + 2)£(s + 2) 

]ys+2 


1 'y 

2 s (s + l)(s + 2)C(s + 2) + s + 


0 ). 


In the region 5fis > —2 the function g is meromorphic with a removable singularity at s = — 1 
and a pole CV = h'( 0 ) at s = 0 , where 

iV s+2 (l + 7 s) 


h(s) = 


2 S (s + l)(s + 2)£(s + 2) 


A direct calculation gives 


C N = h'{ 0) = 


N 2 


- . log N + 7 - log 2 - - - 

2C(2) V 7 g 2 C(2 )J 


We seek to change the contour of integration from a = (To to the contour F consisting 
of the five line segments s = do ± it (t > T), s = o ± \T (—1 < a < uo), s = — 1 + it 
(—T <t<T ), getting 

pcro+ioo 

(4.6) 


^ /*<T0+lOO ^ /* 

— / fif(s) ds = Cat + — / c/(s) ds. 

27T1 d CTn _ ioo 27T1 _/ r 


It remains to show that the contribution of the integral on T is small. Note first that 
|C(s)| = |C( S )I gives \g(s)\ = |g(s)|. As a result only the case 9s > 0 will be considered 
next. Using standard estimates on ( (cf., e.g., :3] > 0 > 13) we have 


(4.7) 
and 

(4.8) 


/ OO 

|fif(cr 0 + it)| 


dt <C 


O' 0 


N 2+a ° log t dt N 2+a ° 

^cro.e 


2 CT °i 3 


f cto rcro ]y2+cro 

i 1 9(<r + iT )l da < J ^ r3 • |C(1 + + iT)| 


J'2-e 


1 


< 0 


r 3 

_/y2+<To 


(T 1 / 2 log T) log 7 T« CT0 „ 


| C (2 + cr + iT)| 

_/y2+o- 0 


dfj 


'O- 0 ^3 l & & ^' cr 0,e J, 5/2-e‘ 

To estimate the contribution of the integrand on the segment — 1 + it (0 < t < T) we follow 
closely the argument from pp. 216-217. The functional equation 

(2ir)‘ 


c(it) = x (it) c(1 - i() 


C(1 + d) 


and the equality 
yield 


C(d) 


C(l + d) 


C(l + d) 

l r (d)| 2 = 

= |x(d)| = 


X(«) = 


2 T(s) cos f 


7r 


f sinh nt 

1 


' t tanh ^ 


2/ 


_ . posh 

t sinh 7rt CUbI1 2 


2vr 


























12 


Employing also tanht < max{t, 1} (t > 0), we get (independently of T > 

The estimates (HU)-(HU) with, say, T = N 2 conclude the proof. 


<C JV. 


□ 


Note also 


Corollary 4.5. g = T_ (log IV + 7 - ||() + O (^) . 

Proof. This is a consequence of Lemma 14.41 and of the well known formulas 

£ ^ = 2^2) + °( iVlo g A 0’ £^ = ^ + 0(logfV). 


a<N 


a<N 


a C(2) 


□ 


Theorem o now follows from (EH, Proposition roi and Lemma roi 


Remark 4.6. Numerical computations show that the error 0(N) given by Lemma 14.41 on 


s n -c n =J2 


y(a)(N-2a) 2 N 2 


2a 2 


a<- 


2C(2) 


lo gJ v + 7 -|-lo g 2-||) 


may not be optimal. Moreover, the graph of Sn ~ Cn exhibits a surprising regularity (cf. 
Figure Q). This was brought to our attention by the referee, who also kindly provided the 
Mathematica notebook. One could hope to improve the theoretical estimate of the error 
by shifting the segment s = — 1 + it further left to the line $is = —1 — 5. The problem 
however is that the argument 2 + s will enter the critical strip 0 < < 1 where lower 

bound estimates for £ are problematic. 



Figure 1 . The plot of Sn — Cn, N < 3000, and respectively N < 20000. 
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Figure 2. The plot of £ a<iV/2 y(a)( ^ 2a) - , N < 10000. 


5. An application to the distribution of quadratic irrationals 


and^T 

qn(x) 


Let x be an irrational number in (0,1) with continued fraction [ao(x), a±(x ),... ] 
be its n th convergent. A classical result of P. Levy states for almost every x that 

lim - log q n (x) = J 
n —»oo n 12 log 2 

The Gauss map T : [0,1) —> [0,1) defined as T(0) = 0 and T(x) = y — [y] if x / 0, has 
the well-known properties 

T([ai, a 2 ,... ]) = [a 2 ,a 3 ,...], 
xT(x)---T n ~ 1 (x) = (-1 ) n (xq n (x) - p n (x)). 

When a; is a quadratic irrational, it is well-known that the limit j3{oj) of ^ log q n (ui) 
exists, and is called the Levy constant of ui. Let AX 2 + BX + C be the minimal integer 
polynomial of u and A = B 2 — 4AC. The length of u> is defined as p(u>) = 2 log £o(w), where 
£q{(jo) = \{uq + wo\/A) is the fundamental solution of the Pell equation u 2 — At; 2 = 4. 

We are interested in the set 7 Z of all purely periodic quadratic irrationals, aiming to 
evaluate 


vr 0 (A)= £ 1 




oo 


uj&TL 

p{u)<X 


Following ^], one defines for each such u = [ai,..., a n ] with n = per(cu) the quantities 

{ n, if n = per(o;) even, 

2 n, if n = per(cu) odd, 

M(lo) = M{a\) ■ ■ ■ M(a n ), 

M( u), if n even, 


m{oj) = n + = 


M(u) 


if n odd. 
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According to (U, Proposition 2.2, we have 

e 0 M = ujT(lo)T 2 (lu) ■ ■ • 
p(uj) = 21oge 0 (w). 

Writing A = / 2 Ao for some fundamental discriminant Ao and some positive integer /, 
one considers the group of units of 0/\ 

u + v\J~A 


A A = 


: (u, v) G Z , u 2 — Av 2 = ±4 


with fundamental unit > 1, in the quadratic field A = Q(\/Ao) endowed with Q-valued 
norm A f and trace tr. One also considers the subgroup 

u + v\fA 


El = 


: (u,v) G Z 2 , u 2 - Av 2 = 4 


of the totally positive units, which is generated by 

ea if A/'(ea) = +1, 
e 2 A if AT{£a) = -!• 


e A = 


In jH] , Section 2, the explicit isomorphism 

A^ ■ F > P'Ai A. 


a b 
c d 


= av + d 


between the hxed point group 

Fu = {g e G = GL 2 (Z)/ ± I : gu = UJ} 
and A a was studied. The inverse of A w acts as 


a: 


f u + vy/~A\ 

u—Bv 

2 

-Cv 

V 2 ) 

Av 

u+Bv 

2 


It sends e A = Qiu> + Qi-\ with i = eper(u;), to fl + = M(uj). Moreover, one has 

AT o \ u = det and tr oA^ = Tr . 

This implies in particular that £o(u) coincides with the spectral radius R(M(co)) of M(u), 
thus 

p(u) = 21ogi?(M(o;)). 

Denote by TZ(A) the set of reduced quadratic irrationals of discriminant A and consider 
the set 

T(N) = Uk,u) : uj £ 1Z = J 11(A), tr (e+(u;) fc ) < n\. 

A>0 2 
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As seen in Section 2, the latter has cardinality T ev (A r ). As shown in f|, Proposition 4.3, 
the map given by 

2m 

eper(u;) 

defines a one-to-one correspondence between Wev(A) and T(N). Denote 


UJ 


where u = [ai,..., a 2m ], 



r(N) = Y 1 = 7r 0 (21oglV). 

s+(u)<N 

Then the identification between W ev (A r ) and T(N) plainly implies as in the proof of jjj, 
Proposition 4.5, the inequalities 

(5.1) Y, r ((N-\y) <* ev (N)< Y 

l<k<2 log N ^ ' l<fc<21og7V 

From the first inequality we infer 


(5.2) 


r(N) < T ev (lV + 1) <C N 2 . 


From (15.IF (15.2F and Proposition II.21 we derive 


<mao < E r(ivi) < E M(* + i)‘) 


l</c<21og N 


l<k<2 log N 


(5.3) 


= T ev (lV+l) + 0 


Y nI 

\ 2<k<2 log N 


= ^ ev (N + l) + 0(NlogN) 
N 2 log 2 


2C(2) 


+ o E (m +£ ). 


The estimate 


r(N) 


N 2 log 2 

2C(2) 


+ o £ (m +£ ) 


follows now immediately from (fT2ll and (15.3F This completes the proof of Proposition ED 
by taking N = e" 2 ". 
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